Effects of a nonlinear perturbation on dynamical tunneling in cold atoms. 
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We perform a numerical analysis of the effects of a nonlinear perturbation on the quantum dynam- 
ics of two models describing non-interacting cold atoms in a standing wave of light with a periodical 
modulated amplitude A(t). One model is the driven pendulum, considered in ref. [1], and the other 
is a variant of the well-known Kicked Rotator Model. In absence of the nonlinear perturbation, 
the system is invariant under some discrete symmetries and quantum dynamical tunnelling between 
symmetric classical islands is found. The presence of nonlinearity destroys tunnelling, breaking the 
symmetries of the system. Finally, further consequences of nonlinearity in the kicked rotator case 
are considered. 



PACS numbers: 05.45.-a, 03.65.-w 

Tunneling is one of the most typical features of quan- 
tum mechanics, concerning oscillations between states 
that cannot be connected in the classical Hamiltonian 
dynamics by real trajectories. The original formulation 
of the problem involved states separated by potential bar- 
riers: the relevant features are semiclassically explained 
in terms of complex solutions of Hamilton equations for 
one-dimensional systems [2] , while a proper treatment in 
higher dimensions is considerably subtler even when inte- 
grability is preserved [3] . Recently a novel kind of tunnel- 
ing, involving transitions between classically separated 
regions in the presence of a non trivial structure of the 
phase space, has attracted considerable attention, both 
theoretically [4-7], and experimentally [8,1,9]. In partic- 
ular a large fraction of these papers focus on physical set- 
tings realized with cold atoms in optical potentials. The 
recent widespread interest and experimental activity in 
Bose-Einstein condensation [10] suggest to check whether 
the presence of Gross-Pitaevskii nonlinearities [11] deeply 
influences the characteristic features of dynamical tun- 
neling. Such a question was already raised and analyzed 
in the framework of the kicked oscillator [12], where it 
was observed how the nonlinear terms typically destroys 
quantum effects induced by symmetry [13,14]. The pa- 
per is organized as follows: we firstly give a few details 
and fix notations for the class of models we are going 
to consider, then analyze two cases: a driven pendulum 
and the kicked rotator; in the last section we briefly re- 
consider pioneering work made on the nonlinear kicked 
rotator [15,16] and supplement it with novel results for 
the quantum resonant case. 



I. GENERAL SETTING. 

We will consider models described by the following 
hamiltonian 

H = Ho + V{>&,i) = ?-+A(t)cos-&, (1) 



where p and $ are canonical momentum and position 
coordinates respectively, expressed in scaled dimension- 
less units, and A(t) is a periodic function of time. Un- 
der appropriate choices of the function A(t), the quan- 
tum version of this model describes an ensemble of non- 
interacting cold atoms in presence of a standing wave 
with a periodically modulated amplitude A(t). The con- 
nection between the scaled variables p and the physical 
ones i9',p' is $ = 2fc^i9 and p — kp jlkjji ( where kr, is 
half of the wave vector of the standing wave and k is the 
effective Planck constant [$,p] = —ik) [17]. 



FIG. 1. Poincare section of an ensemble of classical par- 
ticles for the driven pendulum. The initial ensemble is uni- 
formly distributed in \p\ < 10, 9 £ [— 7T, 7r[; the values of p and 
9 are taken at intervals multiples of the period r. 

We consider two possible choices for the periodic func- 
tion A(t): 



• The first model is the driven pendulum. The ampli- 
tude A(t) is that reported in [1], i.e. — 2a cos 2 (7rf); 
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the period r of time modulation is 1 and the pa- 
rameter a depends on some physical quantities of 
the system, kept fixed in experiments, such as the 
ac Stark shift amplitude, the electric field strength 
and the dipole momentum of the atom. 

• The second model is a variant of the well-known 
Kicked Rotator model [18,17], where the ampli- 
tude A(t) is a sum of 5 functions of period r, 
k cos $^2™^°° S(t — tit). The parameter fc mea- 
sures the kick strength. 
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FIG. 2. (a) The momentum distributions at times t — 
and (b) t — 56r for one rotor with zero quasi-momentum. The 
values of the parameters are a = 10.0, N = 128, L — 500. (c) 
The correspondent tunnelling oscillations of the first moment. 
Time is measured in number of periods. 

Owing to the time dependence of the amplitude A(t), 
in the Kicked Rotator time can be treated as a discrete 
variable, measured in intervals of the period r, and the 
evolution is discrete; in the driven pendulum instead, 
time is a continuous variable. 

The two models present some analogies. In both mod- 
els the potential V($, t) is periodic both in time and space 
with period r and 2?r respectively. Moreover, the two 
models share some symmetries: time-reversal (($,£>, t) — > 
(0, -p, -t)) and parity t) (-0, -p, -i)). 

The quantum Hamiltonian operator of the unper- 
turbed system is obtained from the Hamiltonian (1) re- 
placing the classical canonical variables with the corre- 
spondent operators $ and p — —ikd/di). 

Space and time periodicity are conveniently dealt with 
by using Bloch Floquet Theory. As the potential com- 
mutes with spatial translation of 2n, the quantized mo- 
mentum gets eigenvalues on a discrete lattice p = Tin + /3, 
n G Z and where /3, called the quasimomentum, is the 
analogous of the Bloch phase in solid state physics [19]. 
In the present paper we choose fi — 1 (see the former 



discussion about choice of units), so that the quasimo- 
mentum takes values in the interval [0, 1) (the analogous 
of the first Brillouin zone). 

Quasi-momentum is preserved during quantum evolu- 
tion and parameterizes a single quantum rotor, fixing the 
origin of the discrete lattice in momentum space. The 
evolution of an ensemble of non interacting atoms can be 
modelled by a superposition of the evolutions of indepen- 
dent rotors. 

Because of the invariance of the quantum evolution op- 
erator under time translation by an interval r (U (t+r) = 
U(t)), an evolution operator U T over one period can be 
defined, called Floquet operator. The eigenvalues of 
Floquet operator are e~ 1 ^ , where the real quantities 
e, independent from r, are called quasi-energies. The 
quasi-energy spectrum is invariant under translation over 
2tt?i/t. Quasi-energy plays the same role of energy in 
systems with continuous time variable. 

Concerning the discrete symmetries we mentioned, 
note that changing the sign of p means to change the 
sign of the integer part of p ([p] = n — > —n) and to re- 
place the fractional part {p} of p by 1 — {p} ((3 — > 1 — 0). 
Therefore only rotors with (3 = or (3 = i are invariant 
with respect to these discrete symmetries; rotors with 
(3^0,^ mix each others: a rotor labelled by (3 is mapped 
to another rotor labelled by 0' = 1 — /3. 

Now we discuss some feature of quantum dynamical 
tunneling for the two cases we selected, also considering 
the effect of Gross-Pitacvskii nonlinearities. 
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FIG. 3. (a),(b),(c) Three Floquet eigenstates involved in 
the dynamics of the system; (d) the Fourier spectrum of (p(t)). 
The dashed horizontal line determines the limit of the domi- 
nant frequencies: the five higher peaks emerge. 



II. THE DRIVEN PENDULUM. 

Firstly we analyze the unperturbed system (1) with- 
out the nonlinear term, which is the theoretical setting 
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for experimental data reported in ref. [1]. 

Note that, in the experimental results reported in [1], 
the momentum is measured in units of Ik^h-, which cor- 
respond to the scaled momentum p divided by fc ~ 2.08. 

The classical trajectories are obtained by a numerical 
integration of the Hamilton equations with fourth-fifth 
order Runge-Kutta method. 

In Fig.l the Poincare surface of section for positions 

and momenta p of an ensemble of classical particles is 
shown for the value a = 10.5: the time evolution of 1000 
initial conditions, uniformly distributed in the square 
\p\ < 10, "d G [— 7r, 7r[, is recorded at 100 subsequent mod- 
ulation periods. The symmetries of the system under the 
transformations i9 — > — $ and p — > — p arc clear. A pair 
of stable fixed points of period two (i9 = 0',Po = ±8.65) 
are present. The stability islands, formed by regular tra- 
jectories surrounding the fixed points, are time-reversed 
images of each other. 
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FIG. 4. The momentum distributions at times t — and 
(b) t — lOr in a system of 101 /3-rotors, with f3 equally spaced 
in the interval 0.4 < /3 < 0.6. The values of the parameters 
are a — 10.5, N = 128, L — 500. (c) The correspondent tun- 
nelling oscillations. Time is measured in number of periods. 

Since the Hamiltonian operator is non autonomous, the 
exact Floquet operator of a quantum rotor, with fixed 
quasi-momentum (3, is given by a Dyson expansion. In 
the numerical simulation we use the lowest order split 
method [20]. The Floquet operator is approximated by 
an ordered product of evolution operators on a small in- 
tervals of time At = t/L, with L integer: 

L> T (/3)=Tcxp- t /o T (^+^ t )) 
= rife=i exp _I 4r(™+^) 2 exp 



(2) 



exp 



We use a finite base of dimension N: the discrete 
momentum eigenvalues belong to the finite lattice p = 



(m — y) +/? an d t ne continuous angle variable is approx- 
imated by ■& = t^(to - 1) with m £ Z, 1 < m < N. 

We start by considering the case of a simple rotor, with 
quasimomentum (3 = 0. 

In fig. 2 the evolution of the first moment and the cor- 
respondent momentum distribution are shown for a = 
10. The initial state is a coherent state centered in one 
of the classical symmetric islands: 



f(p) = Ccxp- rpdo cxp- 52( - p - p °f , 



(3) 



where the constant 6 is y/l/2 and C a normalization 
factor. In the simulations we use h — 1, i?o = and 
Po = 8.525. 
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FIG. 5. (a)(p(t)} for u = 5, a = 10, JV = 128; (b) 
(p{t))u — (p(t))o versus t, for u = 3(solid line),7(dashed 
line),10(dotted line) Time is measured in number of periods. 

It can be seen that the maximum of the distribution 
probability oscillates periodically between the two sym- 
metric values +po and — po- In fig. 2(b) the probability 
density is shown at time t = 56r corresponding to a neg- 
ative value of (p): its maximum is peaked at n = —9. As 
we already remarked, for (3 — the quantum system is 
invariant with respect to parity P and time inversion T 
and dynamical tunnelling between the symmetric classi- 
cal stable regions is present, marked by quasi-periodical 
oscillations of the first moment between symmetric pos- 
itive and negative values (approximately (p) w ±6). For 
values of quasimomenta [3 different from or 1/2, exact 
quantum symmetries are broken and periodical oscilla- 
tions are damped and then suppressed. The damping of 
oscillations is faster for values of (3 far from and 1/2. 

The time evolution of the first moment can be ex- 
panded in terms of Floquet eigenstates as: 

<P(*)> = E c i (°K (°) exp- * t( ^ "^E n Xj (n)xt (») 



(4) 
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where the Floquet eigenfunctions \j (p) are expressed in 
momentum representation, the phases (fijj are the cor- 
respondent eigenvalues (connected to quasi-energies by 
<f> = re) and the squared modulus of the coefficients 
c j(0) = (XjlV'(O)) gives the overlap of the eigenfunction 
|Xj) with the initial state. As it can be seen from (4), 
the Fourier frequencies ujj.i, in which the periodic mo- 
tion can be decomposed, are separations between two 
Floquet eigenvalues (ujj.i = cfij — <fin). The dominant fre- 
quencies Ujj correspond to Floquet eigcnstates that more 
contribute to the dynamics, having an high overlap with 
the initial state. 




FIG. 6. The separation time T between P u and Po, 
when \P U (T) — P u= o(T)\ exceeds the value A versus 
u. Starting from above, the six curves correspond to 
A = 0.2, 0.1, 0.05, 0.02, 0.01, 0.005. In the inset a log-log plot 
is shown. 

The Floquet eigenfunctions and eigenphases are eval- 
uated by a numerical diagonalization of the operator. 

By using a finite basis (the momentum is discretized 
using N points), the Floquet operator is reduced to a 
finite matrix N x N and it is calculated as follow. The 
columns of the Floquet operator in the momentum repre- 
sentation are obtained by evolving over one period r the 
momentum eigenstates. Owing to the finiteness of the 
basis, the evolution of the momentum eigenstates ini- 
tially localized on the edges of the basis is affected by 
errors. To overcome this problem, we calculate the Flo- 
quet operator using a basis of double dimension N = 2N, 
so that the variable p takes discrete values in the interval 
p e [-N+1, N] . Then we extract from the matrix N xN 
a nonunitary submatrix of dimension TV in which p varies 
in the range [— ^ + 1+, ^]. The disadvantage of using a 
nonunitary matrix is to find some eigenvalues inside the 
unitary circle. Nevertheless, the number of nonunitary 
eigenphases and the errors can be reduced by increasing 
the dimensions of the bases N and N. 

The pairs of eigenfunctions (Xj>Xi) that dominate the 
dynamics of the system are selected by three conditions 
[21]: a) the maximum presence probability p inside the 
region where the classical stable islands lie (i.e. |n| < 10); 



b) the maximum overlapping probability of each eigen- 
function with the initial state, Pj — (x 3 : |^(0)) | 2 ; c) the 
maximum mutual overlapping probability between the 
two eigenfunctions, PjP%- In fig. 3(a), (b),(c) three Flo- 
quet eigenstates verifying the conditions pj > 0.5, Pj > 
0.07 and P 3 P t > 0.0049 are shown. 

Fourier analysis of (p(t)), with a resolution of 
27r/512 w 0.0123, reveals five dominant frequencies (u)\ = 
0.0859; u> 2 = 0.1841; w 3 = 0.3559; w 4 = 0.4909; w 5 = 
1.1167). A part of the spectrum around zero is repre- 
sented in fig. 3 (d). The dominant frequencies are the 
peaks higher than the dashed horizontal line. 

The differences between the eigenphases of the eigen- 
functions shown in fig. 3 ((f> a = 4.4579, (fib — 4.1035, <fi c = 
3.9198) correspond to the frequencies of the spectrum, 
marked by arrows: u>2 = (fib — <fic, w 3 = <fia — 4>bi 

UA= (fia- <fic- 

Up to now we have considered quantum evolution cor- 
responding to a fixed quasimomcntum (that furthermore 
respects quantum symmetries): we now consider the evo- 
lution of a distribution of rotors with different quasimo- 
menta (3, each evolving with the operator (2) with a fixed 
quasi-momentum (3 [22]: ifi(p,0) = f{p)S n ^. The initial 
momentum distribution f(p) is a coherent state peaked 
in the center of one of the two classical islands of fig.l, 
as in (3). 

During the evolution each /3-rotor evolves indepen- 
dently and the mean value of the observables of the sys- 
tem is an average over different rotors. The time evo- 
lution of the momentum distribution and of the mean 
value of p has been calculated. The data for a = 10.5 
are shown in fig. 4. The spreading in quasi-momentum is 
A/3 = 0.2 around the value (3 = \, which preserves the 
symmetries. The oscillations of (p(t)) correspond to the 
tunnelling of the quantum state between the two sym- 
metric islands in the classical phase space. Note that, 
as in the experimental data [1], during time evolution, 
the momentum distribution of the system maintains its 
maximum localized in the starting island, in the region 
of positive momenta (see the momentum distribution in 
fig. 4(b), corresponding to a minimum value of the aver- 
age momentum). Therefore, in consequence of the aver- 
age over different rotors, the oscillations of (p{t)) do not 
reach negative values, for which higher values of proba- 
bility density in regions of negative momenta are needed. 

Moreover the spreading in quasi-momentum can repro- 
duce the decay of the oscillations in time. In accordance 
with [7], we have also verified that the decay of oscilla- 
tions is faster if the dispersion in f3 is increased. 

As it can be seen from short-period oscillations in 
fig. 4(c), the behaviour of (p(t)} is a superposition of tun- 
nelling oscillations with slightly different frequencies, be- 
cause the spreading in quasimomenta corresponds to a 
spreading of the dominant frequencies of the motion. 

We now study the influence of a Gross-Pitaevskii non- 
linearity in the dynamics: where Schrodinger operator 
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is modified by adding a nonlinear term dependent on 
the squared modulus of the wave function of the sys- 
tem, ii | -0 1 2 . The approximate expression of the evolution 
operator becomes: 



+00 



am 



exp 



exp 



-t^«|Vn| 



(5) 



breaks the sym- 



The multiplicative factor exp" 
metries of the system. 

In presence of the nonlinearity term, the tunnelling os- 
cillations of the first moment loose periodicity in time and 
are progressively suppressed at long times (see fig. 5(a)). 
In fig. 5(b) the difference between the first moment of 
the unperturbed system and that of the perturbed one 
is shown for three values of the nonlinear parameter u: 
u = 3 (solid line), u = 7 (dashed line), u = 10 (dotted 
line). 
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FIG. 7. Classical phase space of the Standard Map (8) for 
K — 2.21 with islands of accelerator modes. The phase space 
is symmetric under space and time inversion. 

We have also calculated the time-averaged first mo- 
ment P{t) — \ Y^=a(p{t')) f° r ni ty values of u equally 
spaced in the interval < u < 10 and we have compared 
the value of P u with Pq in absence of nonlinearity. 

In fig. 6 the time T at which the quantity \P U (T) — 
P u =o(T)\ exceeds the fixed value A versus the nonlinear 
parameter u is shown; each curve corresponds to a dif- 
ferent value of the parameter A from 0.2 to 0.005. The 
intervals of time T decrease algebraically with respect to 
u with the law T = B/u A . In table 1 the power-law ex- 
ponents A and the constants B are shown for different 
parameters A. 



H= ^- + fccosi?^5(t-nr) 



(6) 



n=0 



The correspondent classical map is the well-known 
Standard Map [23], defined on the cylinder: 



Pn+l =Pn + k cos i?„ 
'n+i = #n + T„+i modfa) 



(7) 



After the scaling of the variable p (p = rp/2) and the 
introduction of the parameter K = rfc/2, the map can be 
symmetrized and reduced on the 2-torus [0, 2tt[x [— n, n[: 



Pn+l 
0+, 



where the signs 
the n-th and n - 



t-i ' 

Pn 



■d+ + v 

f K cos j mod(2Tr) 



(8) 



'n+l 



P 



n+l 



mod(2Tr) 



and — refer to instants after and before 
1-th kicks. 



The classical system depends only on the parameter 
K. For K > K c w 0.97 the system undergoes a transi- 
tion to global chaos, the last K.A.M. curve breaks and 
unbounded diffusion in action space takes place. Nev- 
ertheless even for K >> 1 some small stability islands 
survive in the classical phase space, correspondent to ac- 
celerator modes [23-251. 




FIG. 8. Husimi distribution of the state vector at times 
(a)t=800r, (b)t=900r, (c)t=8400r. (d) Tunnelling oscilla- 
tions of (p(t)}. Time is measured in number of periods. 



III. KICKED ROTATOR MODEL. 



The Hamiltonian of the classical system is: 



In fig. 7 the classical phase space of the map (8) for 
K = 2.21 with the stability islands of accelerator modes 
is shown. 
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by h with respect to the classical ones (r = T c {h and 
= k c i/h). Therefore the quantum parameter r plays 
the role of 7i; the classical limit of the system in obtained 
by the limits r — > 0, fc — > oo and K = kr = const. 

For /3 = or /3 = i the system is invariant respect to 
parity P and time inversion T, as in the former case. 
Therefore the Floquct cigenstates belong to invariant 
subspaces with respect to the discrete symmetries; for 
example they can be classified in two classes: odd or 
even with respect to parity [29,30]. 

In the following we fix the value of /3 = and analyze 
first the kicked rotator for a resonant value of t and then 
for a generic value of r. 



5-01 45.05 45.09 
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A. Resonant case. 



FIG. 9. (a) A section of the eigenphases spectrum vs 
k of the Floquet operator. Three eigenvalues for each 
values of k are represented. (b),(c),(d) The momentum 
distribution of the three correspondent eigenfunctions for 
k = 45.02. (e),(f),(g) The three correspondent eigenfunctions 
for k = 45.05. (h),(i),(l) The three correspondent eigenfunc- 
tions for k = 45.08. 

There are two periodic orbits of period two: one formed 
by the pair of fixed points (ir, ±0.7656) and the other by 
the pair (0, ±2.3714). The points of each periodic orbits 
are symmetric respect to p = 0. They are surrounded 
by stability islands inside which the dynamics of the sys- 
tem is regular; each island is delimited by K. A.M. curves, 
which cannot be crossed. During the evolution, ensem- 
bles of points initially located inside the islands centered 
in (7T, ±0.7656) do not mix with those located inside the 
islands centered in (0, ±2.3714). 

The invariance of the system under space reflection and 
time inversion can be seen clearly from the structure of 
the phase space. 

The quantum version of the system is a variant of the 
Kicked Rotator Model [26] , described by the Floquet op- 
erator: 

U T = RKR = e- 4 T(«+/3) 2 g-^cos^^ft+ZJ)^ ( 9 ) 

(the kicked rotator corresponds to quasimomentum ft — 
0). 

For values of r = 4ttM/A (with M, N <E N) the kicked 
rotator undergoes a quantum resonance [27], where the 
spectrum acquires a band structure, yielding ballistic 
transport (see also [28]). 

Having set Ti = 1, the parameters r and k are scaled 



Under the resonance condition r = 4nM/N, the quan- 
tum system is periodic of period N, if M or N is even 
[31]. Therefore the quantum system can be reduced on a 
torus and its Floquet operator becomes a unitary finite 
matrix of dimension N x N. Exact eigenfunctions can be 
calculated. 

The discrete momentum eigenvalues are p = n = in — 
N/2, with m integer, varying in the interval 1 < m < N. 
To make a comparison between the quantum system and 
the classical one on the torus, the classical variable p has 
to be rescaled by the factor 2/t, i.e. p = ^p. Note that 
in the following figures the apex of the quantum variable 
p is omitted. 

In fig. 8 the first moment of p and the Husimi func- 
tion of the momentum distributions at different times are 
shown for M = 1,N = 128 and k = 4.42/r (r w 0.098 
and k = 45.022). For the chosen value of r, the stable pe- 
riodic orbits of the classical map (8) correspond to orbits 
(it, ±15.6) and (0, ±48.3) in the quantum phase-space. 
The chosen initial state is a coherent state (3) centered 
in ($o = 0,Po = - * 0.766 « 15.605). Periodic tunnelling 
oscillations between the two symmetric islands centered 
in ±p take place. 

The tunnelling oscillations are not suppressed even for 
long times. The calculation of (p(t)} in fig. 8(d) is carried 
on for 10 6 modulation periods. 

Note that at t = 8400t w 824.668, when (p(t)) assumes 
a maximum value, the quantum state is mainly localized 
in the island to the upper bound of the torus, centered 
approximately at p w ±48.3098 (see fig. 8(c)). As al- 
ready said, the classical evolution of ensembles of points 
located in the islands centered in p ~ ±2.4(p w ±48.3) 
is independent from the evolution of points inside the is- 
lands centered in p ~ ±0.8 (p w ±15.6), in one of which 
the initial wave packet is localized. 
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FIG. 10. (a),(b),(c) The Husimi distribution of three eigen- 
functions of the Floquet operator, (d) A section of the Fourier 
spectrum of (p(t)). 

The quantum evolution instead couples structures that 
are independent in the classical system. In fact some 
eigenf unctions of the Floquet operator, involved in the 
dynamics, have high probability in both the pairs of is- 
lands (see fig. 10 (a),(b),(c)). 




FIG. 11. Suppression of periodical oscillations of the first 
moment (p(t)} in the presence of nonlinear perturbation ((a) 
u = 5, (b) u = 10, (c) u = 20). The difference A(p(t)) 
is plotted. The values of the parameters are r = 4-7r/128, 
k = 4.42/r, N = 128, L = 2000. 

In fig. 9(a) a portion of the spectrum of the Floquet 
eigenphases versus k is shown. An avoided crossing can 
be seen, close to the value of k used in the calculations 
of this section, i.e. k = — w 45.022. The corresponding 
eigenf unctions were selected by having a support in the 



region \n\ < 25 greater than 0.4. The region \n\ < 25 
corresponds to the classical portion of the phase space 
\p\ = 5= 1-227, in which the stability islands centered 
in pq w ±0.7656 lie. Owing to the invariance properties 
of the system, the eigenfunctions localized in regular re- 
gions of the classical phase space occur in doublets of 
states with opposite symmetries and nearly degenerate 
eigenvalues; wave packets localized in the symmetric is- 
lands of stability are formed by symmetric or antisym- 
metric combinations of such doublets. In fig. 9 (c),(d) an 
example of a quasi-degenerate doublet for k = 45.02 is 
shown; the Husimi function of one state of this doublet 
is plotted in fig. 10(b). 

The dynamical tunnelling can be explained by a three 
states model [30]. In fig.9(b)-(l) the three eigenfunctions 
involved in the tunnelling process for three different val- 
ues of k are shown. Fixing a value of k, the three states 
that take part in the dynamical tunnelling are a doublet 
of quasi- degenerate states with opposite symmetry, lo- 
calized in regular islands, and a third state localized in 
the chaotic region outside the stability islands ((b) for 
k = 45.02, (e) for k = 45.05,(1) for k = 45.08). It is the 
chaotic state that enhance dynamic tunnelling between 
stability islands. 

By varying k, the chaotic eigenstate (fig. 9(b) for k = 
45.02) mixes itself with the state of the doublet sharing 
the same symmetry (fig. 9(d) for k = 45.02) until a com- 
plete exchange between the two states happens (fig. 9(h) 
and fig.9(l) for k = 45.08) [6,30]. 

Beside the triplet of states shown in fig. 9, there are 
other Floquet eigenstates that contribute to the be- 
haviour of the first moment (4) at fixed k; the separations 
between their eigenvalues correspond to the frequencies 
found by the decomposition of (p(t)) in Fourier compo- 
nents. 

The frequencies mostly contributing to the motion 
are revealed by peaks of the power spectrum of (p(t)), 
represented in fig. 10(d) for the parameter k 45.022. 
The spectrum has been calculated with a resolution of 
Auj = 2tt/2 18 ps 0.000024. The three frequencies marked 
by arrows are wi = 0.000057±0.000012, u 2 = 0.000321± 
0.000012, uj 3 = 0.001040±0.000012; u>i, uj 2 and uj 3 corre- 
spond to the difference between eigenphases of three dou- 
blets of quasi-degenerate eigenfunctions: u>i = 4> a — <t>a' = 
5.501869-5.502808,^2 = <?V -<?V = 2.449073-2.448753 
and lu 3 = <\> d - <j) c , = 4.733879 - 4.732833, where 
0a , 4>b, 4>c are the phases of the eigenfunctions plotted in 
fig. 10(a), 10(b), 10(c) respectively. These eigenfunctions 
have a support p in the region |n| < 25 greater than 
0.45 and an overlap probability with the initial state Pj 
greater than 0.01 (P j P l > 0.0001). 

Also in this case we may consider the effect of a non- 
linear perturbation u\ip\ 2 in the Hamiltonian operator. 
The effect of this nonlinear perturbation on the quantum 
dynamics of the Kicked Rotator has already been studied 
in ref. [15,16] but only in non resonant cases. 
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t(.T) 

FIG. 12. Second moment vs time, measured in numbers of 
periods for different values of the nonlinearity: starting from 
above, u = dashed-dotted line, u = 5 dotted line, u — 10 
dashed line, u = 20 full line. The values of the parameters 
are r = 4n, k = 2.5, N = 16384, L = 2000 

Since the perturbation depends continuously on time 
through the wave function ip of the system, the time evo- 
lution operator R between the kicks is approximated as a 
product of evolution operators on small intervals of time 
t/L (L is the number of small steps): 



L 

R= JJe- < * B V i 4*l**l a (10) 

i=l 



The effect of a nonlinear perturbation is to break the 
symmetries of the system thus destroying tunnelling. In 
fig. 11 the difference between the first moment of the un- 
perturbed system (u = 0) and that of the perturbed one 
A{p(t)) = {p(t)) u - (p(*))o is plotted for u = 5, 10, 20. 

We have also analyzed the effect of the nonlinear per- 
turbation on the second moment. In analogy with what 
has been found for the kicked harmonic oscillator un- 
der the resonance condition [14], in the presence of the 
nonlinear perturbation the growth of the second moment 
is slower than in the kicked rotator, even if it remains 
ballistic ((p 2 (i)) = c((3)t 2 ). In fig. 12 the second mo- 
ment vs time is shown: the four curves correspond to 
u = 0, 5, 10, 20 starting from above. This persistence of 
resonant behaviour once the nonlinear term is switched 
on (at least on our observation time scale) is somehow 
surprising, and we hope further theoretical analysis will 
reveal its significance (i.e. if it not suppressed on longer 
time scales). 



B. Generic case. 




01 23456789 



t(«T) 

FIG. 13. Suppression of dynamical localization of the time 
averaged second moment Pi nt (t) in the presence of nonlinear 
perturbation (starting from below, u — 0, full line, u = 5, 
dashed line, it = 10 dotted line, u = 20, dashed-dotted line). 
The parameters values are T = l,fc = 2.5,iV = 128,L = 2000. 

We consider the kicked rotator for r = 1 and analyze 
the effect of the nonlinear perturbation on the first and 
second moment. For this value of r the system is far from 
the classical limit and displays quantum localization. 
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FIG. 14. (a) Tunnelling oscillations in the non resonant 
Kicked Rotator; (b),(c),(d) suppression of dynamical tunnel- 
ing of the first moment (p(t)) in the presence of nonlinear 
perturbation ((b) u = 5, (c)lt = 10, (d)w = 20). The differ- 
ence A{p(t)) is plotted. The parameters values are r = 1, 
k = 4.42, N = 128, L — 2000. 
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FIG. 15. Time T at which |A(p(i))| becomes greater than 
A versus the nonlinear parameter u. Each curve correspond 
to a different values ol A: A = 0.4, 0.2, 0.1, 0.05 starting from 
above. 

In accordance with [16], the nonlinear perturbation 
suppresses localization and give rise to an anomalous dif- 
fusion with an exponent approximately equal to 2/5. In 
fig. 13 a bilogarithmic plot of the time-integrated second 
moment P? nt (t) = \ E^p 2 ^')) for k = 2 - 5 and for dif- 
ferent values of u is shown; the straight line in the figure 
has a slope of 2/5. 

Tunnelling oscillations of (p(t)) for the Kicked Rota- 
tor are plotted in fig. 14(a). In Fig. 14 (b),(c),(d), we 
show the suppression of tunnelling oscillations of the first 
moment for different values of the nonlinear parame- 
ter (b)it = 5, (c)u = 10, (d) u = 20. The differences 
A(p{t)) = (p{t)) u - {p{t)) a are plotted. 

In Fig. 15 the time T at which |A(p(i)}| becomes 
greater that A versus u is plotted. Each curve corre- 
sponds to a different values of A. Note that the evolution 
of (p(t)) in the kicked rotator is calculated at time inter- 
vals equals to t, therefore the estimate of the time T is 
less precise than that calculated for the driven pendulum. 



IV. CONCLUSIONS. 

We have analyzed different physical systems for which 
dynamical tunneling occurs: in particular the driven pen- 
dulum and the kicked rotator. These systems are rele- 
vant for experimental settings recently realized with cold 
atoms. We have stressed several new features: from the 
role of how initially quasimomentum states are assem- 
bled, to the effect of Gross-Pitaevskii nonlinearities, in 
particular revealing subtle features happening in the res- 
onant kicked rotator case. 
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A 


B 


A 


U.UUo 


U.U / lydzU.UUlz 


U.lo4 ( ztu.UUyo 


0.01 


0.0836±0.0012 


0.1649±0.0084 


0.02 


0.0979±0.0011 


0.1664±0.0072 


0.05 


0.1219±0.0011 


0.1908±0.0059 


0.1 


0.1456±0.0010 


0.2031±0.0050 


0.2 


0.1766±0.0012 


0.2193±0.0042 



TABLE I. Fitting parameters of T. 
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